Field theoretic formulation of a mode-coupling equation for colloids 
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The only available quantitative description of the slowing down of the dynamics upon approaching 
the glass transition has been, so far, the mode-coupling theory, developed in the 80's by Gotze and 
collaborators. The standard derivation of this theory does not result from a systematic expansion. 
We present a field theoretic formulation that arrives at very similar mode-coupling equation but 
which is based on a variational principle and on a controlled expansion in a small diniensioneless 
parameter. Our approach applies to such physical systems as colloids interacting via a mildly 
repulsive potential. It can in principle, with moderate efforts, be extended to higher orders and to 
multipoint correlation functions. 



When a suspension of polymer particles is cooled or 
compressed, a rapid slowing down of the dynamics oc- 
curs, and the suspension gradually becomes solid on ex- 
perimental time scales, without any apparent change in 
structure [1]. This colloidal glass transition is reminis- 
cent of the phenomenology of molecular glasses. How- 
ever colloids are conceptually simpler to analyze: the 
interaction potential often has a simple repulsive char- 
acter (instead of a Lennard-Jones form) and their effec- 
tive dynamics is Brownian (instead of Newtonian) . They 
are also experimentally simpler to probe, since colloids 
are much larger than molecules in simple liquids. An 
important class of colloids arc those who interact via a 
bounded, repulsive potential. These, due to the exis- 
tence of a finite energy scale in the potential, exhibit 
a re-entrant behaviour at high density -the glass melts 
upon increasing the density- and the non-interacting liq- 
uid is recovered in the limit of infinite density [2-4]. All 
particles evolve in a thermal bath (the solution) and thus 
undergo individual Brownian motions, while also inter- 
acting via a given pair-potential v. To make our ap- 
proach explicit, we chose to study the harmonic spheres 
model, where the pair-potential v is taken to be of the 
form: v{r) = £ (l — ^) 6* (l — ^), but most of the rea- 
soning will be carried out for an arbitrary, sufficiently 
well-behaved fonction v. This model was introduced by 
Durian [5] in the context of foam mechanics, where v{r) 
plays the role of an effective interaction potential that 
arises from a coarse-graining procedure, but experimen- 
tal realizations in colloids [6, 7] exist, and it became a 
model system to study glassy structure and dynamics [8] . 

The position (t) of each of the A'' particles composing 
the colloidal suspension evolve under Brownian dynam- 
ics, encoded in the following Langevin equations: 



dfj 
dt 



it) = -J2 Vr,v in{t) - r,{t)) + e,W, (1) 



where is a Gaussian white noise with variance 2T (T 
is the bath temperature). It is our goal to obtain quan- 
titative predictions for the dynamics of the dense liquid 



phase of such colloidal suspensions upon approaching the 
glass transition. The only successful first-principles the- 
ory to this day is the Mode-Coupling Theory (MCT) de- 
velopped by Gotze and collaborators [9, 10]. This is a 
closed, self-consistent equation of evolution for the relax- 
ation of density fluctuations in equilibrium supercooled 
liquids which was initially applied to particles evolving 
under Hamiltonian dynamics, but was later extended by 
Szamel and Lowen [11] to interacting Brownian particles. 
No significant difference between these descriptions [12] 
emerges, at least within the MCT approximation. In 
both frameworks MCT predicts a strict dynamical ar- 
rest: below a critical temperature, density fluctuations 
are prevented from relaxing at long times, and ergodicity 
is spontaneously broken. While successes and failures of 
MCT are now well documented [13], a systematic way 
of improving this approximation scheme to overcome the 
listed pitfalls is still lacking, since the original kinetic 
formulation of MCT involves physically motivated, but 
mathematically ill controlled approximations bearing on 
high-order correlation functions, and it contains no a pri- 
ori small parameter. 

The purpose of this letter is to present a new derivation 
of an MCT equation that bypasses several known pit- 
falls at the same time: non-interacting particles (v — 0) 
are exactly dealt with, there exists a small dimensionless 
parameter, the strength of the potential e/T, it follows 
from a well defined variational principle, it can easily be 
extrapolated to higher-orders, and calculations for four- 
point quantities as well as for sheared systems can simply 
be implemented. 

In order to gain insight into what could be a second 
order MCT, the idea of resorting to a field theoretic for- 
mulation is very appealing, since one can then exploit the 
standard toolbox of diagrammatic expansions and ap- 
proximations developped in hard condensed matter and 
particle physics. Several crucial steps have been made 
over the past ten years in this direction. Preliminary 
works [13, 14] have soon been shown to be inconsistent 
with micro-reversibility, a property which is not auto- 
matically conserved by standard approximations in field 
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theory. MCT predicts an ergodic-nonergodic transition 
and one must make sure it does not result from a sym- 
metry breaking approximation. Further attempts [15- 
17] have considerably progressed into the conservation of 
micro-reversibility, but technical difficulties led to either 
a non closed equation for density correlations, or to non 
physical behavior of the solutions to the equations. 

In recent years, Kawasaki and Kim [18] obtained a re- 
sult consistent with reversibility, which led to the same 
equation as that of the original MCT, but this result 
stems from a very cumbersome calculation, giving little 
hope of extending this result to higher orders. In the 
present letter, we suggest to further exploit the many- 
body theory tools used in condensed matter, by formally 
treating our classical particles as bosons. We will see that 
this approach automatically solves several of the prob- 
lems encountered in previous attempts of the formulation 
of a field-theoretic MCT, and provides a transparent way 
to carry the approximations to next order, or to extend 
the calculation to different quantities, such as four-point 
correlators, or to non equilibrium settings, such as in 
sheared systems. We now proceed with a step-by-step 
presentation of our approach. 

The N coupled Langevin equations Eq. (1) can be de- 
scribed by a Fokkcr-Planck equation governing the evo- 
lution of the probability P{{fi},t) of finding each par- 
ticle i at a given position at a time t. As a con- 
sequence of micro-reversibility, the Fokker-Planck equa- 
tion converges towards a Gibbsian equilibrium distribu- 
tion P({rJ,i oo) e-^'<^ "(*'"-*'^'/^. A standard 
result [19] shows that the knowledge of the equilibrium 
distribution allows one to render the Fokker-Planck oper- 
ator Hermitian in the proper basis (this is often called the 
Darboux or supersymmetry transformation). Hermitic- 
ity allows one to interpret this new Fokkcr-Planck equa- 
tion as a quantum mechanical problem for interacting 
bosons. From there we use standard methods of quantum 
field theory [20] to describe the dynamics of the system, 
which is encoded in the following action: 



5 [a, a] = 



a dta + Va • Va 



+ VMaa]/T. (2) 



which is expressed in terms of a pair of complex and 
conjugate fields a and a. The kinetic term reflects the free 
diffusion of particles, and the two-body interactions with 
potential v are expressed, in the quantum formulation, 
by the effective potential that now contains not only two 
but also three-body interactions as follows: 

Vcsip] I p{x)V^v {x - y) p{y) ■ V sv {x - z) p{z) 

Jt,x,yS 



p{x)Asv {x-y)p{y), 



(3) 



t,x,y 



where p = aa is the physical density and v is the pair 
potential between the colloids. 



It is important to notice that the symmetrization corre- 
sponds to a change of basis, so that this field theory does 
not represent directly the physical problem anymore. It 
was shown long ago [21, 22] that the micro-reversibility 
of usual dynamical field-theories obtained from the reg- 
ular Fokker-Planck equation is represented by a compli- 
cated, non linear transformation, making it very difficult 
to preserve when performing mode-coupling approxima- 
tions [16]. In the symmetrized theory, micro-reversibility 
is simply encoded in the hermiticity of the symmetrized 
operator, which is a symmetry easy to check and con- 
serve even when performing approximations. Further- 
more, setting the pair-potential to cancels the effective 
potential in Eq. (3) and one recovers, without approxima- 
tion, the free diffusion of colloids. Even if the dynamics 
described by the action Eq. (2) is only related to the real 
dynamics by a change of basis, careful analysis of the the- 
ory shows that far from initial and final conditions (in the 
"bulk" of the time window), the difference between the 
modified dynamics and the real dynamics vanishes. Fi- 
nally, keeping in mind that the pair potential v{r) has an 
energy scale e, we sec that this approach gives a satisfac- 
tory basis for a perturbation expansion in powers of the 
dimensionless parameter e/T. Our approach yields a the- 
ory that is expressed with ladder operators a and a that 
do not directly represent the physical density p = aa. 
Introducing it by hand via a Lagrange multiplier, a field 
A, we arrive at a field theory involving four independent 
fields, that we group into a single vector 4> = {o-,a,X,p). 
To conclude the layout for the diagrammatic expansion 
to come, we prefer working with fields defined by devia- 
tions around the saddle of the action, which describe a 
homogeneous and isotropic liquid state of mean density 
Pq. Thus we set = {^,^,0,pq) + (*,*,A,(5p), 
and obtain a four field field theory characterized by the 
following propagator: 
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(4) 



where M(fc) = 
teraction part 



of the 



(1 -f pov{k)/Ty - 1 . The in- 
action is made of two cubic 
terms, one is J' ^{ki,k2,k^)5p(ki)5p(k2)5p(kz) and the 
other is j* X{ki)^{k2)'^{k^). where 7(^1,^2,^3) = 
■:^[ki ■ k2 v{ki)v{k2) + perms], and the symbol J* 
stands for 6{ki + k2 + k^). We now turn to the 

procedure allowing us to determine the correlations of 
our fields. 

The physical quantity that we are ultimately interested 
in is the matrix element of the two point correlator of 
the theory G that describes density-density correlations. 
Since our theory contains four fields, G is a 4 x 4 ma- 
trix, with 10 independent entries. Our goal is to obtain a 
closed equation bearing on the Sp-Sp element only. Per- 
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forming a double Legendre transform of the dynamical 
partition function, one obtains a functional of the corre- 
lator G, called the 2PI or the Luttinger-Ward functional, 
$ [G] . A careful analysis of this functional, which can be 
found in modern field theory textbooks [23, 24], shows 
that it has several remarkable properties. It not only 
provides a variational principle to obtain the correlator 
G (the functional attains its maximum when evaluated at 
the true correlator) but also gives access to the inverse 
of the correlator, since it is obtained as the functional 
derivative of $ [G] with respect to G. Diagrammatically, 
it is composed of all two particle irreducible diagrams 
(2PI), allowing for simple truncations of the complete ex- 
pression of $. Finally, any truncation of $ can be shown 
to preserve the symmetries of the action, which wc use 
to conserve micro-reversibility when performing approx- 
imations. To obtain a self-consistent approximation for 
the two point correlator, one constructs an approxima- 
tion for the 2PI functional by selecting a certain sub- 
class of diagrams that contribute to it. For example, the 
two simplest diagrams that contribute have the following 
topology: 



(5) 



Then an expression for the vertex function E (the inverse 
of the correlator) is obtained by functionally differentiat- 
ing with respect to G. Finally, a self consistent equation 
is obtained by exploiting the relationship that exists be- 
tween S and G (sometimes referred to as the Schwinger- 
Dyson equation): (Gq ^ — I][G]) G = 1. This variational 
approach can be seen as the dynamical counterpart to 
the density functional theory of liquids. 

So far, all these considerations are in principle exact. 
Wc now present the simplest self-consistent approxima- 
tion that can be obtained within this formalism, and we 
will see that we obtain an equation that has the exact 
same structure as that of the Mode-Coupling equation. 
We now exploit that our theory contains a small param- 
eter, namely the strength of the potential, to select the 
lowest order beyond mean-field. We have two vertices in 
the theory, one of order 4 in e, and one of order 0. We 
want to stay to the lowest non-trivial order, so we can 
neglect the former. 

Wc only retain the simplest diagram in the expression 
of the 2PI functional Eq. (5), and will justify this a pos- 
teriori. Wc obtain an expression for the vertex function 
that can be inserted in the Schwinger-Dyson equation to 
yield: 



Go^G{k, t)= / r(fc, q)G{k -q,T- t')G{q, T - t')G{k, t') 



the mode-coupling equation, in which the memory kernel 
is a quadratic functional of the correlators, except that 
it applies to a matrix instead of a scalar. In order to 
write down an equation that involves the density-density 
correlator only, we must express all other correlators in 

terms of G(fc, t — i') = (6p{—k,t')5p{k,t)^. At the mean- 
field level, all correlators are proportional; we use these 
proportionality relations and insert them into Eq. (6). 
The proportionality coefficients involve various powers 
of e/T, and one then verifies that, when inserting the 
proportionality relations into the expression of the 2PI 
functional, all contributions coming from the diagrams 
that we neglected are indeed of higher order in e/T. 

In Eq.(5) only the topology of the diagrams is repre- 
sented, but one has to draw all possible diagrams from 
the vertices of the theory. Even for the simplest water- 
melon diagram, this involves 11 independent diagrams. 
Fortunately, to lowest order (order 2 in e/T) only one 
diagram survives, and the final evolution equation for G 
is: 



= -dlG{k, t) + n{kfG{k, t) 



(7) 



+±. f Mik, q)G{q, T - t')G{k -q,T- t')d,.G{k, t'), 
where the memory kernel has the following expression: 



M{k, q) = <fc{q) + {k~ q)\{k - q) 



(8) 



Note that, as usual in field theoretic formulations, one 
has had to resort to a further approximation, i.e. set- 
ting —v{k)/T = c{k). This results from our treating 
the statics and the dynamics on equal footing. At this 
order of approximation, this replacement is correct, as 
can be seen with a perturbation analysis of the equi- 
librium liquid, and allows for direct comparison with 
the regular MCT result. The statics of equation (8) is 
closely similar to the original mode-coupling equation, 
apart from the slightly different wave-vector dependence 
that the original mode-coupling approach predicts, in 
which the factors q^ and (fc — q)^ in the rhs of (8) are 
replaced with k ■ q and k ■ {k — q), respectively. Assum- 
ing that the density density correlation function does not 
decay to zero at large times, one makes the usual ansatz: 
limt_j.oo G(A:, t) = pQS{k)f{k), where S{k) is the static 
structure factor related to the direct correlation function 
by S{k) = 1/(1 — pQc{k)), and seek an equation for the 
non-ergodicity parameter f{k). By Laplace transform 
methods one easily obtains: 



f{k) poS{k) 



(6) 1 - /(fc) 87r2fc4 



This is a matrix equation in which all correlators appear. 
Note that Eq. (6) in itself already has the structure of 



M{k,q)S{q}Sik - q)f{q)f{k - q) 

(9) 

We then numerically solve this equation with an iterative 
procedure. The only input is c{k) for the equilibrium 
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liquid, than can be calculated e.g. within the Hyper- 
Netted Chain approximation. Exactly as in the case of 
standard MCT, one finds that there exists a transition 
line Ta/ctIp) above which f{k) = is the only solution, 
whereas below Tmct{p), a nonzero f{k) is found where 
ergodicity is spontaneously broken. A posteriori, we are 
inclined to view the MCT equation as a high-temperature 
expansion. In Fig. 1 we show the resulting f{k) at pack- 
ing fraction 0.53 and temperature 10~^. The qualitative 




FIG. 1: Non-ergodicity parameter f{k) calculated with the 
regular MCT kernel (full line) and within our modified kernel 
(dotted line), at packing fraction 0.53 and temperature 10^'*. 

behavior is very similar, except for the limit fc — > 0, where 
Eq. (9)) is found to give f{k) ^ 1. The modification 

fe-i-O 

of the kernel implies that the absolute value of the tran- 
sition temperature is slightly modified. For example at 
the packing fraction 0.53, the regular MCT transitionis 
located approximately at Tmct ~ 10^'*, whereas our 
rough numerical estimate is Tmct ~ 8.10""'. 

In this letter we have presented a comprehensive ap- 
proach to write down mode-coupling equations based on 
a variational principle. In the example of bounded in- 
teractions, we have shown that when the strength of the 
interaction is taken as an expansion parameter, it is possi- 
ble to write down, to lowest-order, a mode-coupling equa- 
tion similar to the regular MCT equation. Our strategy 
can be extended in a variety of directions. The most obvi- 
ous one is retaining higher orders in the expansion param- 
eter e/T. The resulting equation for the 7-dimensional 
order parameter G will pick up a contribution to its 
memory kernel. Retaining, after appropriate substitu- 
tions based on the leading order expansion (6), the next 
order in e/T seems a tedious yet quite accessible task. 
It would also be of interest to examine whether qualita- 
tive differences show up if the full set of ten equations 
(6) were solved. On our to-do list we also have more 
pressing wishes like implementing the so-called "thermo- 
dynamic of histories" formalism [25, 26] and probing the 
relationships between ergodicity breaking and dynamic 
phase transitions. It would be interesting to investigate 



sheared systems and compare our approach with existing 
extensions of MCT [27]. The fate of ergodicity breaking 
and dynamic phase transitions under shear also belongs 
to our open questions. 
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